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To reach the target of smaller pressure drop and better heat transfer performance, packed beds with small 


tube-to-particle diameter ratio (D/d)<10) have now been considered in many areas. Fluid-to-wall heat transfer 


coefficient is an important factor determining the performance of this type of beds. In this work, local flu- 


id-to-wall heat transfer characteristic in packed beds was studied by Computational Fluid Dynamics (CFD) at 


different Reynolds number for D/d,=1.5, 3.0 and 5.6. The results show that the fluid-to-wall heat transfer coeffi- 


cient is oscillating along the bed with small tube-to-particle diameter ratio. Moreover, this phenomenon was ex- 


plained by field synergy principle in detail. Two arrangement structures of particles in packed beds were recom- 


mended based on the synergy characteristic between flow and temperature fields. This study provides a new local 


understanding of fluid-to-wall heat transfer in packed beds with small tube-to-particle diameter ratio. 


Keywords: Packed bed, Fluid-to-wall heat transfer, Field synergy principle, Computational fluid dynamics, Heat 


transfer intensification 


Introduction 


Packed beds have been extensively applied to chemi- 
cal and industrial processes such as reactors, separators, 
dryers, filters and heat exchangers. To reach the target of 
smaller pressure drop and better heat transfer perfor- 
mance, packed beds are now being developed into small 
tube-to-particle diameter ratio (D/d,<10) (Dixon et al. 
[1]). Flow and heat transfer in this type of packed bed 
play an important role in determining the performance of 
the bed, and have therefore attracted an ever increasing 
attention of numerous investigations. A good local qua- 
litative understanding and description of flow and heat 
transfer are important for the design and operation of 
such type of packed bed. 
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For heat transfer in packed beds, studies on fluid-to- 
wall heat transfer in packed beds could be classified as 
theoretical study, experimental research and numerical 
simulation. In theoretical studies, packed beds are gener- 
ally treated as one-dimensional pseudo-homogeneous me- 
dia, two-dimensional pseudo-homogeneous media and 
heterogeneous media. Then fluid-to-wall heat transfer 
coefficient is obtained as a part of governing equation or 
a boundary condition. Hsu and Cheng [2] and Ozgumus 
et al. [3] had made a perfect summarization on two- 
dimensional pseudo-homogeneous model in packed beds. 
The one-dimensional pseudo-homogeneous model now 
could be easily solved out. And, the analytical solutions 
for two-dimensional pseudo-homogeneous model and 
heterogeneous model were also presented by Dixon and 
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Nomenclature 
A; sectional area of an element (m°) 
dy diameter of particle (m) 
D diameter of the packed bed (m) 
h local fluid-to-wall heat transfer coefficient 
: (W/(m? K)) 
A p/L Pressure drop per unit length (Pa/m) 
Pr Prandtl number 
qwlz local heat transfer in the wall (W/m°) 
r radial position (m) 
Re, Reynolds number based on superficial velocity 
Re, Reynolds number based on interstitial velocity 
Ty avelz local average temperature of fluid (K) 
Ta Wall temperature (K) 


Cresswell [4]. However, the local information about flu- 
id-to-wall heat transfer is all absent in theoretical studies. 

Experimental studies on fluid-to-wall heat transfer in 
packed beds mostly focus on developing an explicit cor- 
relation of fluid-to-wall heat transfer coefficient to pre- 
dict its value. Several correlations for fluid-to-wall heat 
transfer coefficient in packed beds can be found in the 
literatures. Both Miroliaei et al. [5] and Jorge et al. [6] 
had made a summarization on correlations for fluid-to- 
wall heat transfer coefficient in packed beds. Wen and 
Ding [7] also summarized correlations of fluid-to-wall 
heat transfer coefficient in packed beds and concluded 
that the Li-Finlayson correlation agrees best with their 
experimental data. However, the effect of local structure 
on fluid-to-wall heat transfer is not easy to be obtained in 
experimental studies. 

Computational Fluid Dynamics (CFD) is extensively 
used to investigate the fluid-to-wall heat transfer due to 
the fact that it is an effective tool for gaining local details 
in packed beds. Nijemeisland and Dixon [8] simulated 
fluid-to-wall heat transfer in packed beds by CFD, the 
results of which are in good agreement with the correla- 
tions proposed in the literatures. Guardo et al. [9] inves- 
tigated the influence of turbulence model in CFD model- 
ing on wall-to-fluid heat transfer in packed beds. Log- 
tenberg and Dixon [10] also used computational fluid 
dynamics to obtain the values of the dimensionless wall 
heat transfer coefficient. Unfortunately, to the best of our 
knowledge, how the local packing and flow structure are 
related to fluid-to-wall heat transfer in packed beds and 
what packing structures are optimal for heat transfer en- 
hancement are still questions in the literatures. 
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AT Temperature gradient (K/m) 
fluid velocity in the entrance of packed beds 
U -1 
(ms) 
y volume of the macro-cylinder with the length 
a of dz 
y volume of particles in the macro-cylinder 
Pd with the length of dz 
x position in x direction (m) 
y position in y direction (m) 
4 axial position (m) 
Greek Symbols 
E average void fraction through the bed 
E void fraction in the infinitesimal cylinder 
c between z and z+dz 
0, synergy angle in a mesh element in the x-y 
! plane 
Dave average synergy angle in the x-y plane 


To address the above issues, the aim of this work is to 
establish a relation between local structure and fluid-to- 
wall heat transfer in packed beds with small tube-to- 
particle diameter ratio. We first use the Discrete Element 
Method (DEM) to obtain three-dimensional packed beds 
randomly packed with spherical particles in a cylindrical 
tube. Then the flow and local fluid-to-wall heat transfer 
in the beds are simulated by CFD. Finally, field synergy 
principle is used to analyze the effect of local packing 
and flow structure on fluid-to-wall heat transfer in pack- 
ed beds, and the optimized structures are given. 


Numerical methodology 


Governing equations 


For the simulations in this study, full three-dimensional 
governing equations were applied. The conservation of 
mass, Momentum and energy are needed for the flow and 
heat transfer of single phase fluid. 

The continuity equation for conservation of mass is 
given as: 


—(pu;)=0 (1) 


where p is the mass density, x; are the Cartesian coor- 
dinates, and u; are the velocity components. 
The momentum equation for conservation of momen- 
tum is given as: 
——(pu,u;)=-P +44 pg, (2) 
Ox ; Ox; Ox; 
where p is the static pressure, g; is the gravitational 
body force, t; is the stress tensor, defined as follows for 
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Newtonian fluid: 


ty =| iu Ou; 2 ou, 
l Ox; = OX; 
where u is the dynamic viscosity, and 6; is the Kro- 


necker delta. 
The energy equation for conservation of energy is 


2 Ou; 
u 


y 3 
oa @) 


given as: 
fa] fa] or ô 
—(u;h)= k; + UT; 4 
a Out) zeua] a n) (4) 
where A is the specific enthalpy, and fk; is the thermal 
conductivity. 


For turbulent flow, as Guardo et al. [9], one-equation 
model is chosen for simulations on turbulent flow and its 
equation is as follows: 


J j J 
(5) 
ĉuj Ou; Oe Kv? 
C 
ee. | Ge, ej | PPT 


where ø, is Prandtl number of turbulence pulsation 
kinetic-energy, cp is empirical coefficient, and yi, is tur- 
bulent pulsating viscosity, defined as follows: 


ec pel (6) 
where c’, is empirical coefficient, and / is the turbulent 
pulsation length. 


Geometric Model 


The packing structure has an important influence on 
the flow and heat transfer characteristics in packed beds. 
Therefore, the establishment of it should approach to the 
real situation as far as possible. This study developed a 
three-dimensional packed bed structure with uniform 
spherical particles randomly packing in a cylindrical tube 
using Hertze-Mindlin non-linear contact model in Dis- 
crete Element Method (Ucgul et al. [11]). The packed 
beds with D/d, of 1.5, 3.0 and 5.6 are shown in Fig. 1. In 
the simulations, the diameter of particles is 3mm and the 
length of the cylindrical tube is about five times the di- 
ameter of particles. Although the length of the beds is 
short, the entrance effect only affects magnitude of heat 
transfer and not impacts the morphological characteris- 
tics of heat transfer. The material of spherical particles is 
glass and that of tube wall is steel. As to coordinate sys- 
tem of the bed, the z direction is parallel to the axis of the 
bed, with the x and y directions determined by right-hand 
screw rule. For the bed with D/d,=5.6, local void fraction 
of DEM simulation were compared with MRI data of 
Sederman et al.[12] and predicted values of correlations 
of Klerk [13] and Mueller [14], as shown in Fig. 2. It can 
be seen that the distribution of void fraction along the 
radial direction of the bed from simulation are, as a 


whole, in agreement with MRI data and predicted values 
of correlations, though there are a little discrepancies 
between simulation and data/correlation, especially at the 
maxima and minima. Local void fraction of the simula- 
tion is slightly greater than experimental data because the 
simulated bed is loose packing compared with the expe- 
rimental packing. Even so, the packing structure built 
using DEM could reflect the real situation of packed beds, 
considering the deviation of experiment and simulation 
and the change of real height of packed particles in the 
radial direction of the bed. 


(a) Didy=1.5 


(b) D/d,=3.0 (c) Did,=5.6 


Fig. 1 Geometry models of packed beds with tube-to-particle 
diameter ratio of 1.5, 3.0 and 5.6. 
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Fig. 2 Comparisons of radial distribution of void fraction 
between CFD calculation and MRI measurement and 
prediction of correlations in the packed bed with 
D/d,=5.6. 


Mesh generation 


Before mesh generation, the zones near contact points 
of particle-particle and particle-wall must be modified for 
the good mesh quality and the calculating convergence in 
numerical simulations. In this work, particles were con- 
tracted by 0.5% of their diameter. As to this method, Di- 
xon et al. [1] and Nijemeisland and Dixon [8] had made a 
conclusion that both the 99.5% and 99% sphere size 
models show negligible difference from the touching 
model in aspect of velocity distribution. Triangular mesh 
was chosen for the surface of particles and tube wall and 
tetrahedral mesh for fluid zone, as shown in Fig. 3. The 
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suitability of mesh used for the simulation is checked and 
the packed bed with D/d,=1.5 is selected for mesh-inde- 
pendence verification. Three sets of meshes with total 
element number of 1683408, 3217115 and 7132803 are 
used for the test and the maximal mesh sizes in the fluid 
zone are of d,/10, d,/20 and d,/40 correspondingly. Fur- 
thermore, the maximal mesh size on the particle surfaces 
is limited to d,/30. It is shown in Fig. 4 that the deviation 
of velocity magnitude in meshes with the maximum 
mesh sizes of d,/10 in fluid zone is very small with that 
of other two meshes. The mesh with the maximum size 
of d,/10 in fluid zone is chosen for all packed beds, in 
order to avoid the big computer hardware requirements 
and shorten the calculation time. 


K SEH Seep 


(a) Shell mesh 


(b) Volume mesh 
Fig.3 Meshes in the packed bed with D/d,=1.5. 
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Fig. 4 Variations of velocity magnitude with mesh numbers in 
the packed bed with D/d,=1.5. 


Simulation Methodology 


Simulations on flow and heat transfer in packed beds 
were carried out in ANSYS FLUENT 12.0. Three-dimen- 
sional Navier-Stokes equations, namely Eqs.(1—2), were 
firstly used to simulate the laminar flow of fluid through 
the packed bed. The one—equation turbulence equation, 
namely Eq.(5), was added for the turbulent flow. In this 
study, three interstitial Reynolds number Re, of 10, 100 
and 1000 were chosen to simulate the flow and heat trans- 
fer in packed beds. According to the transition regime of 
Dybbs and Edwards [15], the flow is laminar for Re,= 10 
and 100, and turbulent for Re,= 1000. Air was chosen as 
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the fluid through the bed and its physical property parame- 
ters are: p = 1.225 kg/m’, c, = 1006 J/(kg:K), ky = 0.0242 
W/(m: K) and u = 1.78 x 10° Pa-s. The operation condi- 
tion pressure is 101325Pa. Boundary conditions were set 
as velocity-inlet for inlet of packed bed, pressure-outlet 
for outlet of packed bed, no-slip wall condition for sur- 
face of particles and wall. The pressure-velocity coupling 
was carried out by the SIMPLE scheme. Second-order 
upwind schemes were used for the convective terms in 
the momentum equation. Simulations were not stopped 
until the residuals fall below 1 x 10° for density, velocity 
and turbulent kinetic energy. Moreover, the drag in the 
surface of particles was as a monitor for simulations of 
flow. Then the predicted flow results were used to simu- 
late fluid-to-wall heat transfer by solving the energy equ- 
ation, namely Eq.(3). Boundary conditions were set as 
373K for inlet of packed bed, heat insulation for surface 
of particles, and 293.15K for tube wall. Second-order 
upwind schemes were used for the convective terms in 
the energy equation. Simulations were not stopped until 
the residuals fall below 1 x 10° for energy. The average 
temperature in the surface of particles was also as a mon- 
itor for simulations of heat transfer. 


Validation 


To validate the reliability of CFD simulation on packed 
beds, the CFD results of pressure drop were compared 
with the predicted values by correlations of Ergun [16], 
Zhavoronkov et al. [17] and Reichelt [18], as shown in 
Fig. 5. It can be seen that the simulated results are in 
good agreement with those estimated using Zhavoronkov 
et al.’s and Reichelt’s correlations for Re,= 10 and 100. 
The pressure drop in packed beds with small D/d, is un- 
derestimated by Ergun’s equation at low Reynolds num- 
ber due to wall friction or viscose forces effects at the 
wall. The simulated results are slightly greater than those 
estimated by correlations for Re, = 1000, which may 
caused by the entrance effect and global shrinking of 
particles. What also can be seen from Fig. 5 is that pres- 
sure drop per unit length increases with the increasing of 


Re,=1000 


- eW- - Simulations —o— Reichett [18] 
103 —4—Zhavoronkov et al. [17] _—*— Ergun [16] 
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Fig. 5 Comparison of values of pressure drop between CFD 
calculation and prediction of correlations. 
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tube-to-particle diameter ratio due to the decreasing of 
void fraction. However, the increasing rate of it decreases 
with the increasing of tube-to-particle diameter ratio pos- 
sibly on the account that the decreasing rate of void frac- 
tion also does. 


Results and Discussion 


Local fluid-to-wall heat transfer coefficient 


Local fluid-to-wall heat transfer coefficient could be 
obtained through dividing temperature difference be- 
tween the average temperature of fluid and wall temper- 
ature in infinitesimal cylinder by the heat rate transferred 
through the wall between z and z+dz. The calculating 
equation could be expressed as follows: 


= - Ww = 7 (7) 
f ave dy w 
Fig. 6 shows the distributions of fluid-to-wall heat 
transfer coefficient and void fraction along the bed at 
different Reynolds number for various D/d,. Local void 
fraction ¢, in infinitesimal cylinder between z and z+dz 
could be obtained using the following equation: 


V 
e, = 2 (8) 
Va 


vA 


h 


Zz 


What can be seen from Fig. 6 are as follows: i) fluid- 
to-wall heat transfer coefficient is periodic oscillating 
along the bed with wave crest and trough alternately pre- 
senting. The fluid-to-wall heat transfer coefficient attains 
to the maximum at the inlet of the bed, then quickly de- 
creases, and finally becomes periodic oscillating. 

The distribution of fluid-to-wall heat transfer coeffi- 
cient along the bed is similar to that of local void fraction, 
just with their wave crest and trough in different posi- 
tions. In packed beds with same D/d,, local fluid-to-wall 
heat transfer coefficient increases with the increasing of 
interstice Reynolds number, as well as difference value 
between peak and trough adjacent. And, the positions of 
wave trough in the distributing curve move to the outlet 
direction of the bed with the increasing of Re,. Under the 
condition of equivalent interstice Reynolds number Re,, 
local fluid-to-wall heat transfer coefficients decrease with 
the increasing of tube-to-particle diameter ratio D/d,, as 
well as different values between peak and trough adja- 
cent, which may be dependent on the decrease of oscilla- 
tion in void fraction near the wall. The three rules above 
will be explained in detail using field synergy principle in 
the next section. 

Fig. 7 shows temperature contour maps in the x-z 
plane for different D/d, at Re,=100. It can be seen that 
periodic thermal boundary layers present near the wall, 
which further confirms the existence of the oscillation of 
wall heat transfer coefficient along the bed. In addition, it 
also can be found that the temperature gradient in the 
region away from the wall is much less than that near the 
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Fig. 6 Local fluid-to-wall heat transfer coefficients along the 
bed at different Reynolds number and local void frac- 
tion along the bed 


wall for the cases studied. It elucidates that the mixing 
heat resistance in the region away from the wall is much 
less than the heat resistance near the wall for the cases 
studied. In other words, the periodically oscillating beha- 
vior of fluid-to-wall heat transfer coefficient along the 
bed is, to a great extent, from the periodical change of 
void fraction near the wall. The reason why the oscillat- 
ing characteristic of heat transfer coefficient along the 
bed cannot be obtained in theoretical models may be that 
void fraction is set to constant in the model, but void 
fraction, especially near the wall, is oscillating along the 
bed with small tube-to-particle diameter ratio. 
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Fig. 7 Temperature contour maps in the x-z plane for different tube-to-particle diameter ratio at Re,=100. 


Field synergy analysis 


Field synergy principle, proposed by Guo and Li [19], 
elaborates the physical mechanism of convection heat 
transfer in the view of the synergy between flow field 
and temperature field. This theory demonstrates that the 
performance of convection heat transfer depends not only 
on the velocity of fluid, the physical properties and the 
temperature difference between fluid and wall, but also 
on the synergy angle of velocity field and temperature 
field. In other words, the better the synergy between ve- 
locity and temperature fields is, the higher is the intensity 
of heat transfer under the same boundary conditions of 
velocity and temperature. 

It can be concluded from the previous section that 
there are some relationship between the oscillation of 
fluid-to-wall heat transfer coefficient along the bed and 
the distributions of local void fraction and local packing 
structure. However, how void fraction and packing struc- 
ture affect fluid-to-wall heat transfer still needs to be 
further accounted for by field synergy principle. In nu- 
merical simulations, the synergy angle 6; in a mesh ele- 
ment between velocity and temperature gradient can be 
calculated from Guo and Huang [20]: 


0; = arcos es (10) 
‘olvz| 

The local average synergy angle avez in the x—y plane 

in the position of z, can be acquired from Guo and Huang, 


[20]: 
(11) 


Fig. 8 shows the distributions of fluid-to-wall heat 
transfer coefficient h, and local average synergy angle 
Oave z along the bed. It can be seen from Fig. 8 that syner- 
gy angle at the inlet attains to the minimum, and thus the 
synergistic characteristics there is the best throughout the 


bed leading to the strongest intensity of heat transfer. 
Then the synergy presents a periodically oscillating cha- 
racteristic, so does the corresponding heat transfer coef- 
ficient. This may be caused by the periodic variation of 
local void fraction near the wall. 

Fig. 9 shows the distributions of velocity and temper- 
ature in the x-y plane for D/d,=1.5 at different Reynolds 
number. It also can be seen that the synergistic characte- 
ristic at the inlet is the best throughout the bed, and the 
synergy as well as the arrangement structure of particles 
presents a periodicity along the bed. Typically, a triangu- 
lar structure formed by three particles and wall (Zone 1), 
periodically presents along the bed, as shown in Fig. 9. In 
this triangular structure, the synergy is good at the corner 
between two particles close to the entrance of the bed, 
whereas a bit poor at other place. It is consistent with the 
distribution of local average synergy angle along the bed 
in Fig. 8. 

The other phenomenon, observed from Fig. 9, is that 
eddy appears if interstice Reynolds number attains to 100, 
which is in favor of the synergy enhancement. The heat 
transfer intensity is directly proportional to flow velocity. 
Local heat transfer coefficients therefore increase with 
the increasing of flow velocity. Moreover, the oscillating 
amplitude of the distributing curve of local heat transfer 
coefficients increases with the increasing of flow velocity. 
This is because the heat transfer is significantly enhanced 
in the location of small synergy angle with the increasing 
of flow velocity, whereas that in the location of large 
synergy angle changes little. 

The synergy characteristics in packed beds with D/d, = 
3.0 and D/d, = 5.6 has the same distributions as in the 
packed bed with D/d, = 1.5, as well as the influence me- 
chanism of flow velocity on heat transfer, as shown in 
Figs.10—13. Just the local average synergy angle increas- 
es with the increasing of tube-to-particle diameter ratio at 
the equivalent interstice Reynolds number. 

The decrease of synergy characteristics with the in- 
creasing of tube-to-particle diameter ratio may be due to 
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the fact that the amount of the triangular structure formed 
by three particles and the wall like Zone 1 in Figs. 9, 11 
and 13 reduces near the wall with the increase of D/d,. 
This causes the decrease of values and axial oscillating 
amplitude of local heat transfer coefficients. Based on 
this rule, we can predict that the oscillating behavior of 
heat transfer coefficient along the bed will disappear only 
if tube-to-particle diameter ratio is large to an extent 
since the uniformity of void fraction prevails near the 
wall in the axial direction of the bed. In Figs. 11 and 13, 
the triangular structure can be observed in a little place 
near the wall, and eddy also appears if flow velocity 
reaches a specific value. In fact, the synergy between 
velocity field and temperature field formed in this trian- 
gular structure, on the whole, is the best among all struc- 


Temperature 


Fig.9 Streamline and temperature contours in the x—z plane in 
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Fig. 8 Local fluid-to-wall heat transfer coefficients and aver- 
age synergy angle along the bed for D/d,=1.5. 
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Fig. 11 Streamline and temperature contours in the x-z plane 
in the packed bed with D/d,=3.0. 


tures near the wall. Therefore, the triangle structure is 
suggested to be arranged near the wall as far as possible. 

Based on the reinforcement mechanism of heat trans- 
fer, the synergy between velocity field and temperature 
field is good at the place where the disturbed flow forms 
(Yang and Tao [21]). In a randomly packed bed, the dis- 
turbed flow forming in two kinds of structures is strong 
in the region away from the wall: one is the rhombus 
structure with a pair of edges having a gap and another 
pair of edges not doing, like Zone 2 in Figs. 11 and 13, 
where the continuously disturbed flow may form; another 
is the rhombus structure with both two pairs of edges 
having a gap, like Zone 3 in Fig. 13, where the cross dis- 
turbed flow may form. The mixing of fluid is enhanced, 
and a jet may produce if the fluid flows through these 
two structures. Based on the existing theory, the cross 
disturbed flow is the best way to strengthen heat transfer 
among these two kinds of disturbed flows (Yang and Tao 
[21]). Therefore, the rhombus structure with both two 
pairs of edges having a gap is suggested to be arranged in 
the region away from the wall as far as possible. 

To validate the reliability of the study on heat transfer 
in packed beds, the average wall Nusselt number ob- 
tained using Eq.(9) in the simulations are compared with 
predicted values of correlations suitable for the cases 
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Fig. 12 Local fluid-to-wall heat transfer coefficients and av- 
erage synergy angle along the bed for D/d,=5.6. 


studied, as shown in Fig. 14. It can be seen that CFD data 
are in good agreement with predicted values of correla- 
tions shown in Tab. 2 with deviation within 25%. 


d, pL 
Nu, = | hy ade (9) 
kL% 


Conclusions 


To understand fluid-to-wall heat transfer coefficient in 
packed beds, in this study, we firstly built randomly pac- 
king structure of particles in packed beds with small tube- 
to-particle diameter ratio (D/d,<10) using DEM, and then 
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Fig. 13 Streamline and temperature contours in the x-z plane in the packed bed with D/d,=5.6. 
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Fig. 14 Comparisons of average wall heat transfer coefficients 
between the simulations and the correlations. 


simulated the fluid-to-wall heat transfer by CFD. Finally, 
we demonstrated the characteristics of local fluid-to-wall 
heat transfer and its relationship with local structure and 
flow using field synergy principle. The following major 
conclusions can be drawn. 

(1) In packed beds with small D/d,, the fluid-to-wall 
heat transfer coefficient is periodically oscillating along 


the bed with wave crest and trough alternately presenting. 


It attains to the maximum at the inlet of the bed, then 
quickly decreases, and finally becomes a periodical os- 
cillating. The distributions of fluid-to-wall heat transfer 
coefficient along the bed are similar to that of local void 
fraction, just with their wave crest in different positions; 

(2) In packed beds with same D/d,, the fluid-to-wall 
heat transfer coefficient increases with the increasing of 
interstice Reynolds number Re,, as well as difference 
values between peak and trough adjacent; 

(3) Under the equivalent interstice Reynolds number, 
the fluid-to-wall heat transfer coefficient decreases with 
the increasing of tube-to-particle diameter ratio D/d,, as 
well as difference values between peak and trough adja- 
cent. Based on this rule, we predict that the oscillating 
behavior of the local heat transfer coefficient along the 


bed will disappear only if tube-to-particle diameter ratio 
is large to an extent since the uniformity of void fraction 
near the wall prevails along the beds; 

(4) The distribution of fluid-to-wall heat transfer coef- 
ficient along the bed and the changing rules of it with 
flow velocity and tube-to-particle diameter ratio are elu- 
cidated through field synergy analysis on velocity and 
temperature fields. The triangular structure formed by 
three particles and the wall is suggested to be arranged 
near the wall and the rhombus structure with both two 
pairs of edges having a gap is suggested to be arranged in 
the region away from the wall. 


Acknowledgements 


This research is supported by the National Natural 
Science Foundation of China (51276181, 51476173) and 
the National Basic Research Program of China (2011CB 
710705). 


References 


[1] Dixon, A.G., Nijemeisland, M., Stitt, E.H., (2006), 
Packed tubular reactor modeling and catalyst design us- 
ing computational fluid dynamics, Advances in Chemical 
Engineering, Vol. 248, pp. 307—389. 

[2] Hsu, C.T., Cheng, P., (1990), Thermal dispersion in a 
porous medium, International Journal of Heat and Mass 
Transfer, Vol. 33, pp. 1587—1597. 

[3] Ozgumus, T., Mobedi, M., Ozkol, U., Nakayama, A., 
(2013), Thermal dispersion in porous media-A review on 
the experimental studies for packed beds, Applied Me- 
chanics Reviews, Vol. 65, pp. 031001-1-19. 

[4] Dixon, A.G., Cresswell, D.L., (1979), Theoretical predic- 
tion of effective heat transfer parameters in packed beds. 
AIChE Journal, Vol. 25, pp. 663—676. 

[5] Miroliaei, A.R., Sharhraki, F., Atashi H., (2011), Compu- 
tational fluid dynamics simulations of pressure drop and 
heat transfer in fixed bed reactor with spherical particles, 


170 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


Korean Journal Chemical Engineering, Vol. 28, pp. 
1474-1479. 

Jorge, L.M.M., Jorge, R.M.M., Giudici R., (2010), Expe- 
rimental and numerical investigation of dynamic heat 
transfer parameters in packed bed. Heat and Mass Trans- 
fer, Vol. 46, pp. 1355-1365. 

Wen, D.S., Ding, Y.L., (2006), Heat transfer of gas flow 
through a packed bed. Chemical Engineering Science, 
Vol. 61, pp. 3532-3542. 

Nijemeisland, M., Dixon, A.G., (2001), Comparison of 
CFD simulations to experiment for convective heat 
transfer in a gas—solid fixed bed, Chemical Engineering 
Journal, Vol. 82, pp. 231-246. 

Guardo, A., Coussirat, M., Recasens, F., Larrayoz, M.A., 
Escaler, X., (2005), Influence of the turbulence model in 
CFD modeling of wall-to-fluid heat transfer in packed beds, 
Chemical Engineering Science, Vol. 60, pp. 1733- 1742. 
Logtenberg S.A. Dixon A.G., (1998), Computational fluid 
dynamics studies of fixed bed heat transfer, Chemical 
Engineering Process, Vol. 37, pp. 7—21. 

Ucgul, M., Fielke, J.M., Saunders, C., (2014), Three- 
dimensional discrete element modeling of tillage: deter- 
mination of a suitable contact model and parameters for a 
cohesionless soil, Biosystems Engineering, Vol. 121, pp. 
105-117. 

Sederman, A.J., Alexander, P., Gladden, L.F., (2001), 
Structure of packed beds probed by Magnetic Resonance 
Imaging, Powder Technology, Vol. 117, pp. 255-269. 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


J. Therm. Sci., Vol.25, No.2, 2016 


Klerk, A., (2003), Voidage variation in packed beds at 
small column to particle diameter ratio, AIChE Journal, 
Vol. 49, pp. 2022-2029. 

Mueller, G.E., (1992), Radial void fraction distributions 
in randomly packed fixed beds of uniformly sized sphe- 
res in cylindrical containers, Powder Technology, Vol. 72, 
pp. 269-275. 

Dybbs A., Edwards R.V., (1984), Fundaments of trans- 
port phenomena in porous media, Martinus Nijhoff, Dor- 
drecht. 

Ergun, S., (1952), Fluid flow through packed columns, 
Chemical Engineering Progress, Vol. 48, pp. 89-94. 
Zhavoronkov, N.M., Aerov, M.E., & Umnik, N.N., (1949), 
Hydraulic resistance and density of packing of a granular 
bed. Journal of Physical Chemistry, Vol. 23, pp. 342-361. 
Reichelt, W., (1972), Zur Berechnung des Druckverlustes 
einphasig durchstromter Kugel- und Zylinderschuttungen. 
Chemie-Ingenieur-Technik, Vol. 44, pp. 1068-1071. 

Yagi S., Wakao N., (1959), Heat transfer from wall to 
fluid in packed beds, AIChE Journal, Vol. 5, 1959, pp. 
79-85. 

Li C., Finlayson B.A., (1977), Heat-transfer in packed- 
beds-re-evaluation, Chemical Engineering Science, Vol. 
32, pp. 1055-1066. 

Demirel Y., Sharma R.N., Al-Ali H.H., (2000), On the 
effective heat transfer parameters in a packed bed, Inter- 
national Journal of Heat and Mass Transfer, Vol. 43, pp. 
327-332. 


